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SYMBOLS 


depolarization factor 
k = index where i ® ice and w * water 


A. = - abc 
J 2 


«D 

/ 


ds 


(a^ + S)^''^(b^ + S)^''^(c^ + + S) 


where a, b and c are ellipsoid axres and S is distance along 

the respective axis 

molecular radius, constant 

form factor 

frequency 

relaxation frequency 

Boltzman's constant 

radius of i = ice, a = air, w = water 

where 


w 


r 


a 



0.916 
w 

100 - 

w 




1/3 


temperature 

volumetric fraction of i = ice, w = water 
percent liquid water by weight 
attenuation coefficient 
penetration depth 

dielectric constant of i = ice, w = water, ds * dry snow, 
ws « wet snow, 0 = static value, » = optical limit 
real part of c 
imaginary part of e 
viscosity i 


• extinction coefficient 
e 

■ absorption coefficient 

K » scattering coefficient 

» wavelength in free space 

p = density of s ■ snow, i = ice 

k 

T « relaxation time 
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DIELECTRIC PROPERTIES OF SNOW 


W. H. Stiles and F, T. Ulaby 
Remote Sensing Laboratory 
University of Kansas Center for Research, Inc. 
Lawrence, Kansas 66045 


ABSTRACT 

This paper is a review of the dielectric properties of snow in the 
radio frequency range from 100 KHz to 35 GHz. Applicable dielectric 
mixing formulas are discussed and compared to available experimental 
datu. Below 1 GHz, the dielectric behavior of snow is well understood and 
has led to development of sensors to measure the liquid water content of 
snow, while above 1 GHz, the scarcity of experimental data and the 
discrepancies between the few available measurements, particularly when 
liquid water is present, are responsible for the rather sketchy understanding 
of the dielectric behavior of snow in the microwave region. 


V 


1.0 INTRODUCTION 

As the 1iT9>ortance of snow as a renewable resource becomes Increas- 
ingly recognized, methods of assessing snow properties In a timely fashion 
over large areas become indispensable. It Is clear that, for both timely 
observations or large ai«a1 coverage, the physical handling of the snow 
must be minimized. Remote sensing techniques satisfy the above require- 
ment; however. Interpretation of the remotely sensed data Is not always 
straightforward. These data are usually In the form of processed quanti- 
ties which are determined by the parameters of the measurement system 
and the characteristics of the observed scene (In this case, snow). In 
the radio region of the spectrum, the Interaction of electromagnetic en- 
ergy with snow Is governed by Its geometrical and electrical (dielectric) 
properties; therefore, a more complete description of the dielectric 
properties of snow should lead to a better understanding of remote sens- 
ing observations of snow. 

Several review articles of the electrical properties of ice and snow 
have been published over the last twenty years [1-4], The most recent of 
these reviews, by Glon and Paren [3], stressed the electrical behavior 
in the spectral range that Is particularly suited to glacial sounding 
(DC i;o 100 KHz). This paper presents a review of the dielectric proper- 
ties of snow, with special eir^hasis on its behavior at the radio frequen- 
cies above 100 KHz. Included In this review are: (1) a discussion of 
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mixing formulas used to characterize the dielectric properties of a snow 
mixture In terms of Its constituent parts— air, Ice and liquid water (If 
present), and (2) a summary of experimental data reported In the litera- 
ture. 


2.0 PROPERTIES OF CONSTITUENTS 

The electromagnetic propagation properties of a material are (tefined 
In terms of Its magnetic permeability, u, and Its relative complex di- 
electric constant, c. For most naturally occurring materials. Including 
snow, u • vq. the magnetic permeability of free-space. Hence, the pro- 
pagation Is governed solely by e, which has the complex form: 

c » e‘ - je” (1) 

Since snow Is a heterogeneous mixture of air, ice and, under certain 

conditions, liquid water, the constituent dielectric properties must be 
examined first In order to understand the dielectric properties of the 
mixture. For practical pirposes, the characteristics of air are Indis- 
tinguishable from those of free space. Water and Ice, however, exhibit 
dielect! k behavior which can be described, at least to the first order, by 
the Debye equation: 

tn • t 

e - (2) 

1 + j2lffT 

or In terms of Its real and imaginary parts: 


e’ » e + X 

* 1 + (2TTfx)^ 

£. « 

1 + {2irfT)^ 


(2a) 

(2b) 


2 



where 


cq * f ^ 0 ^ static dielectric constant, dlmsenslonless 

c, * ^ optical limit of dielectric constant, 

dimensionless 

T « relaxation time of the material, s 
f * electromagnetic frequency, Hz 

This equation describes the contribution to the polarizability of a 
polar molecule (H 2 O) from the permanent dipole moment of that molecule. 

An applied electric field tends to align the dipole against the thermal 
forces which Induce disorder or "relaxation" and are a function of the 
temperature and viscocity of the dielectric material. The relaxation time 
as derived by Debye for spherical polar molecules In a viscous medium Is a 
good approximation for water [5]: 

T . iULii (3) 

kT 

where n Is the viscocity, a Is the molecular radius, k Is Boltznisn's 
constant and T Is the temperature. 

Empirical expressions based on experimental results have been de- 
veloped for the constants In equation (2) for water as functions of 
physical temperature [5-8] and are reviewed In [9]. Figure 1 illustrates 
the behavior o^ e' and e" of water at 0"C as a function of frequency. 

The peak In e"^ at 8.5 GHz results from the Debye dispersion being centered 
In the microwave region. The frequency at which the peak occurs 1$ termed 
tte relaxation frequency, f^, where f^ • (2irT)"^. 

For the frequency Interval of emphasis, 100 kHz to 35 GHz, the 
dispersion causes the dielectric properties of water to vary significantly. 
Below approximately 1 GHz, low frequency approximations to Equations 2a 
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•nd 2b yield «'„ ■ cqw " *«w^* ^ 

exact form must be used. For liquids or solids, interactions between 
molecules and viscosity variations tend to violate some of the assumptions 
of the Debye equation giving rise to a non-unique relaxation time. The 
Cole - Cole [10] modification to the Debye form (Equation 2), i.e. , 
raising the J2vfT term to the power (1 - a), allows this equation to 
describe the dispersion for a distribution of relaxation times. For 
water, a has been determined to be less than 0.03 [8] indicating that 
the Debye form for a single relaxation time is adequate, but for ice 
and snow this is not always the case. 

The dielectric properties of ice, at least over a large frequency 
interval, can also be described by the Debye equation of the same form 
as for water. Equation 3. While the relaxation frequency of water lies 
in the microwave region, for ice (as a result of the higher viscosity compared 
to that of water) the relaxation f*^uency is reduced to the kilohertz 
region and as in the case for water, it is temperature dependent. Auty and 
Cole [11] measured the dielectric coefficients for pure ice and found 
that ■ 3.10, that had a minimun value at O'* corresponding to 
f^ ■ 7.5 KHz, and that had a minimum value of 91.5 at -1** C. 

A convenient form for exming the spectral behavior of a material 
describable by the Debye form is the Cole-Cole diagram which is a plot of 
the imaginery part c" versus the real part e'. Curves defined by Equa- 
tion 2 are semi-circles in this representation. Figure 2 from Auty and 
Cole's data [13] illustrates measurements of pure ice and ice with impuri- 
ties along with their respective Debye curves. Each symbol corresponds 
to a measurement at frequencies denoted in KHz. Two observations are 
apparent from this figure: (1) the DC conductivity of sample (b) causes 
the behavior to deviate from the Debye form below 1.5 KHz and (2) since 
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the optical Unit 1$ the left-most Intersection of the semi-circle with 
the e'-ax1s, this representation is not very Illuminating for the high 
frequency response because of the compression of the variation In e". 

It should be noted that In the frequency range below 10 KHz, several 
loss mechanisms are present [3,12,13], primarily due to Impurities, structural 
effects, and molecular or Ionic defects. A good overview of these dis- 
persions Is given In the review by Glen and Paren [3] and In a paper by 
Paren and Glen [12]. These dispersions may be describable by the Debye 
form but do not occur at the same relaxation frequencies or have the same 
strengths. In the high-frequency limit however, these values will each 
approach their optical limits and consequently will have very little 
contribution above 100 KHz. Therefore, the macroscopic dielectric behavior 
cari be described by effective quantities that Include contributions from all 
of the dispersion mechanisms. 

For frequencies much greater than t1^ relaxation frequency In Ice, 
f/f© ■ 2irft^»l, the following high frequency approximations can be made 
to the rationalized form of the Debye Equation (2a and 2b); 


'l * ®1- 


(4a) 


e 


N 

1 


^10 • ^ 1 » 
2irfT^ 




(4b) 


From measurements for which tl^se approximations are valid, between 10 lOlz 
and 1000 ^z, the optical limit has beeri found to be Insensitive to 
teaq>erature (below 0**C) and Independent of frequency. Values from several 
experimental Investigations [1, 14 ?1] for are given in Table 1 and 
excellent agreement generally Is observed; hence, the value of ** 3.15 
from Gumming [14] will be used as being representative. While the real 
part Is well-behaved, the imaginary part Is dependent on both frequeiKy 


and temperature as expected from Equations (3 and 4b). Fif/ure 3 is a 
survey of the available data on as a function of frequency for oure- 
or fresh-water Ice at two temperatures. For this frequency Int^r^al, 
Equation 4b Indicates that should decrease as f“S therefore on a log- 
log scale, the behavior should be linear with a negative slope. The ex- 
perimental data given in Figure 3 shows, however, that the Debye form Is 
Inapplicable above 1 GHz, Indicating that the single resonance assumption 
Is not valid. Evans [1] attributed the observed spectral behavior of 
e”. to contributions by infrared dispersions. 

Since the theoretical expressions for e". are not adequate in the 
microwave frequency range, experimental results must be relied upon even 
though considerable varability exists in the data above 10 GHz. This 
variation is at least in part a consequence of the problems associated 
with making accurate measurements of the loss factor of a very low-loss 
material, which Ice Is. 

3.0 MIXING FORMULAS 

The dielectric properties of a heterogeneous mixture such as snow 
must be described by mixing formulas that incorporate the dielectric 
properties and proportions of the component materials. In the rigorous 
formulations, structure will also be a factor. Snow can be treated as 
either a two-comoonent or three-component mixture consisting of Ice, air 
ana liquid water. The two-component formulation is usually applied to 
dry snow (air and Ice) and sometimes to wet snow through an additional 
appllcati'^n of the formula for dry snow and water. In the three- component 
formulations, air. Ice and water are each accounted for as Individual 
C(xi^ponents. 
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Both theoretical and purely empirical mixing formulas have been 
reported In the literature. In this section, the discussion will be 
limited to theoretical approoches, while discussion of the regression 
models will be presented In the section on dielectric properties. 

The developments of mixing formulas have usually been restricted to 
modeling one component as regularly shaped inclusions uniformly distributed 
within a host or background material. The mathematical difficulties of 
handling Inclusion shapes and Inter-Inclusion interactions have 
required the choice of ellipsoidal or spherical particles. This limitation 
poses no real problem for snow that has undergone melt-freeze 
metamorphism as the particle shapes tend toward spherical anyway, however, 
very fresh snow is not usually even close to spherical. Applicability 
for fresh snow therefore may be suspect but also may be difficult to examine 
because handling requirements and the action of methamorphism immediately 
start modifying its structure towards more stable thermodynamic shapes. 

Mixing formulas which have been applied to snow are given in Table 2. 
With the exceptions of Formulas 2 and 3, which are variations on the 
formulation due to Weiner [1,22] and which attempt to use a single factor to 
describe a three-dimensional shape (and are therefore Inherently unsound), 
all the other formulas can be related to one another. 

Formulas b, 6 and 7 are descriptive of ellipsoidal inclusions with 
axes a, b and c and associated depolarization factors where 
k denotes the inclusion type, either ice or water, and j denotes the axis. 
Usually two of the axes are asswned equal, thus since the sum of the 
Aj must equal one, only one Integration is needed. Formula 5 is the 
most rigorous In that (the effective dielectric constant in the 
Immediate region of the particle) can be different from either of the 
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constituent!; or the mixture average. This formulation Is the most general; 
however* It Is also the most complicated to apply. Formula 7 makes the 
approximation that the dielectric constant of the mixture, 

and appears to provide an adequate fit to the available data. Both of these 
formulations assume that the dry snow dielectric properties are already 
known, while Formula 6 treats the three component case. In applying Formula 
6 to measured dielectric data. Col beck [24] found that the ellipsoidal Ice 
particles may be approximated as spherical In shape without incurring 
significant error in the results, but for’ the water Inclusions the shape 
factor was found to be very Important. 

Formula 1, which applies to dry snow only. Is a special case of 
Formula 5; if the particles are spherical (i.e., = Ag * A^ = 1/3), 

Formula 5 reduces to Formula 1 upon replacing and with and , 
respectively, and setting Similarly, Formula 4 can be derived 

fran Formula 5 for dry snow conditions. 

The expressions given by Formula 9 is special case of a more general 
formulation derived by Tinga et al. [26] for ellipsoids of dielectric 
covered by a layer of material of dielectric e 2 contained in a host material 
of dielectric ^s given In Table 2, Formula 9 is for spherical Ice 
prticles covered by a uniform layer of water. Formula 8 is the same as 
Formula 9 with no water present (I.e., dry snow). While Colbeck [24] states 
that the water distribution in the snowpack is not in spherical layers around 
the ice particles. Formula 9 has been found to given acceptable results when 
compared to experimental data, as discussed in a later section. 

4.0 DIELECTRIC PROPERTIES OF SNOW 

The dielectric properties of dry and wet snow are discussed In 
separate sections because the addition of liquid water to dry snow sig- 
nificantly alters its properties. 
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4.1 Dry Snow 

Because the relaxation frequency of ice is in the KHz range, dry 
snow is not expected to exhibit a dispersion behavior above 100 KHz. 

Although surface effects and other processes may cause additional low- 
frequency dispersions [5], the mixture with air, which can be considered 
to have f ^ “ makes the relaxation frequency of snow always higher than 

that of ice. The relaxation frequency for natural dry snow remains below 
approximately 100 KHz. Dielectric measurements, conducted for den- 

_3 

sities ranging between 0.08 and 0.41 g an and for temperature between 
-10®C and 0°C, show that the relaxation frequency of dry snow, in 

the 5 KHz to 60 KHz range [1,12,28]. 

Ambach and Denoth [25] measured the frequency dependence of and 
the loss tangent, tan 6^^, for fine-grained (<0.25 mm radius) and coarse- 
grained (>0.25 mm radius) dry snow. The loss tangent is defined as the 
ratio of the imaginary part to the real part of the dielectric constant. 

Their results, given in Figure 4, indicate different relaxation spectra 
for the two samples; the coarse-grained sample exhibits a single Debye 
type relaxation, while the fine-grained sample exhibits a wide relaxation 
spectrum indicative of a range of relaxation frequencies. This result does not 
imply that the size is the effect; the depolarization factors of the ice 
particles were very much different form the two cases: fine grained A<0.03 
and course grained A>0.20. Above 10 MHz, the real part is independent of 
frequency and the effect of grain-shape on the relaxation spectrum in the 
high frequency limit is negligible. 

For frequencies above 1 MHz, the dominant physical parameter that 
determines the dielectric constant of dry snow is its density. An alternate 
and useful expression for density is the volume fraction of ice , which 

is related to the ice density and the snow density through; 
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( 5 ) 


Figure 5, from Gumming [14], depicts the dependence of on at 
9.375 GHz, along with curves defined by mixing Formula 4 (In Table 2) 
and several empirical expressions [3,25,29]. Although this data set 
was measured at 9.375 GHz, It Is applicable above 1 MHz where Is 
frequency- Independent. Cimmlng [14] also found to be Independent 
of temperature as expected from Ice and air properties. The empirical 
fits [3,25 and 29] are approximately what one would get from a volunetricly 
weighted average of the dielectric constants of the constituents. Any of 
the mixing formulas applicable to dry snow will also provide a good 
fit to this data. One useful simple formulation for the real part results 
from Formula 4 with = 3.15, 

e'ds = (1 + 0.508 (6) 

While Is Independent of both temperature and frequency, c''^^ 

Is sensitive to temperature as shown In Figure 6 and to frequency by 
Inference from Figure 3 for Ice. Although spectral measurements of 
dielectric properties have been made up to around 100 Wiz, In the micro- 
wave region single-frequency measurement programs have predominated and precise 
determination of the spectral dependence is not available. As a result, 
one must rely on mixing formulas, which often lead to Inconsistencies. For 
example, Royer [2] found that If Formula 2 of Table 1 Is used, different 
values of F are required to fit the real and Imaginary parts of to 
experimental data suggesting that a better model Is needed and therefore 
Formula 2 should not be used. 

An expression derived from Formula 8 (of Table 1) by Ulaby et al. 

[9] Is given by 
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( 7 ) 


0.34 V.e% 

J.H , 1 I 

(1 - 0.417V^)^ 

and provides a fit close to Cunning's [14] data. 

For dry snow. If one avoids the low frequency dispersion region, 

the real part of Is well-known and easily determined from density and 
the mixing formulas; for the Imaginary part, however, the fits are 

approximate due to the limited accuracy of the available data on the Imaginary 
part of the dielectric constant of Ice. 

4,2 Wet Snow 

As 0^*0 1s approached, liquid water appears In snow and causes a sig- 
nificant change in the dielectric properties of the mixture becaure both 
the real and imaginary parts of the dielectric constant of water are 
much greater than those of ice. 

10 MHz - 1 GHz Region . For the frequency interval between approxi- 
mately 10 MHz and 1 GHz, propagation loss through dry or wet snow is very 
small. Consequently, most of the research activities conducted in this 
frequency range have concentrated on the dependence of g' . on V , the 

Ww W 

volume fraction of water. The use of mixing formulas is facilitated by 
the fact that for f > 10 MHz, e'^ = 3.15 and, as mentioned earlier, it 
is frequency and temperature independent, and for f < 1 GHz, e' = e'wn* 
the static dielectric constant of water, which, of course, is frequency 
independent and weakly dependent on temperature. Hence, the dielectric 
constant of the ice-air-water mixture is also frequency-independent and 
approximately temperature-independent. 

The empirical expressions given In Figure 5 for dry snow have been 
extended to wet snow by adding a tenn that accounts for the water volume 
fraction, V^, 

'"ws * ’ * * •’''w (8) 
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where the constants (a,b) were found experimentally to have the values 
(2.22, 23.6) by Howorka [In 3] and similar values (2.205, 21.3) were 
determined by Ambach and Denoth [25]. 

The above expression provides a useful approximation of the 
relationship between e' and V , although a more exact representation 

Infd W 

using Formula 7 of Table 1 allows for modeling the water particles as 
ellipsoids. With the assumption that the depolarization factors 
and A 2 are equal, implying that A^ = 1 - 2A^, Ambach and Denoth [25] 
employed experimental data to compute A^, by fitting the data to Formula 
7. Their results are shown in Figure 7. It is observed that A^ is a 
function of V^; however, Aj s 0.05 seems to be a good average depolarization 
factor. 

I detailed modeling of liquid distribution, Colbeck [24] 

suggested that water actually appears in two basic shapes (fillets and 
veins) because of the geometry of wet snow ice-grains. Snow that has 
undergone melt-freeze cycles tends to form multiple crystal clusters and 
the fillets and veins of water form in the air spaces between three or 
two grains, respectively. 

The practical application of dielectric measurements in this fre- 
quency-range is in the measurement of liquid water, which will be covered 
in a later section. 

Above 1 GHz . In the microwave region, both real and imaginary 
parts of are frequency-dependent, therefore frequency extrapola- 
tion questionable. Four sets of data have been reported in the 
literature for e as a function of snow liquid water content V . 

Wb Inf 

Gumming [14] reported the variation in loss-tangent in response to V^, 
however, his data covered a narrow range extending between » 0 and 0.006. 
Figure 8 presents dielectric measurements given by Sweeney and Colbeck 
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[23] at 6 GKz and Tobarlas, et a1. [30] at 9.4 GHz, together with the 
results reported by Linlor [31] for both frequencies. Also shown are 
curves computed using the mixing formulas Indicated on the figure. The 
curves due to Linlor [31] were computed from expressions he had developed 
on the basis of experimental data acquired over the 4-12 GHz region. 

These expressions are: 


4 

e' = 1 + 2p + bV 
ws ^s w 

(9) 

« 

II _ 1 .0994a r —\ — 
ws f ws 

(10) 


where 

b = 5.87 X 10'^ - 3.10 x 10‘^f - 4)^ 

a = V [0.045(f - 4) + 0.066(1 + a)] dB/cm (11) 

w 

with f in GHz and with a representing the a\irage grain diameter in mm. 
The ranges of validity are stated to be: 0.1 mm < a ^1.0 mm and 0 < 

< 0.12 [31]. 

The experimental data in Figure 8 show significant differences at 
both frequencies, especially for the imaginary part. In addition, the 
curves generated from mixing Formulas 5, 6 and 9 also show some 
discrepancies. The behavior at 6 GHz shows that Formulas 5 and 
6 give similar results as expected because of their common 
basis. Formula 6 was applied using spherical ice particles (A^.j = 1/3), 
ellipsoidal water particles (A^^ = 0.04, an average fran Figure 7), and 
the constituent dielectric properties. The results from tl Is formula 
fall near an upper bound on the data indicating that at the higher values 
of V^, a different value of A is needed. Formula 9 shows good agreement 
with e' data from Linlor [31] up to about V * 0.1; however, the 
curvature does not follow the correct trend and therefore should only be 
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used with great caution. prediction for e" Is even further In 

W w 

error. 

At 9.4 GHz, again the data from Uni or [31] are much lower than 
the other data. Formula 9 gives a better fit to Linlor's data than at 
6 GHz, while Formula 6 fits the data of Tobarlas, et al. [30] fairly well 
with the same depolarization factors used at 6 GHz. 

The differences In dielectric constants for these sets of data suggest 
that the dielectric behavior of wet snow In the microwave rerlon Is 
Influenced by some parameter(s) besides snow density. Dielectric mixing , 

theory and observations of Ambach and Denoth [36] Indicate that snow grain 
size or structure should not affect the properties. It is the authors’ 
belief that grain structure does Indirectly influence the dielectric 
properties. The method of preparing the samples of reference [23] was to 
flood the snow sample and then drain It. Colbeck [37] states that this 
leads to multi-grain clusters of three and more particles and the liquid 
Inclusions exist at grain boundaries. This results in ellipsoidal liquid 
Inclusions which may explain the better fits obtained by Formulas 5 and 6. 

Colbeck [37] also states that In seasonal snowpack two grain clusters are more 

probable, indicating that the water will take on a ring-like shape. The 

data obtained by Linlor [31] was for snow that had been subjected to "vigorous" 

manual mixing which probably resulted in single ice grains coated with an ^ 

approximately uniform layer of water. Such a geometry lends itself to a good 

fit by Formula 9. Since this is an artificial snow condition, it should be ♦ 

considered at most lower limit for e’ and e" . 

Current understanding of the microwave dielectric behavior of snow 
may be simmarizcd as follows: 

(a) For dry snow, is Independent of frequency, approximately 
temperature- independent, and its variation with snow density Is fairly well 
established (Figure 5) and may be approximated by Equation 7 or Formula 6. 
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(2) For dry snow, e"jg varies with temperature, frequency and den- 
sity. Very few experimental data are available, particularly above 

10 GHz. 

(3) For wet snow, the presence of water adds another dimension to 
the complexity of the dielectric behavior of snow, and yet available 
data on e' - and e‘* . are even fewer than those reported for dry snow, 
thereby resulting in a rather sketchy understanding of the variation of 
e with frequency, and water content. However, Formulas 5 and 6 
appear to provide the best fit to non-mi xed snow. 

Clearly, a vigorous experimental program designed to establish 
the dielectric behavior of snow "'s called for, particularly when free 
water is present. 

5.0 PENETRATION DEPTH 

The penetration depth in a medium is defined as the depth at which the 
average power of a wave traveling downward in the medium is equal to 1/e 
of the power at a point just beneath the surface of the medium. For a 
medium with uniform extinction coefficient ic^, the penetration depth 6p 
is given by 

«p = • 

In the general case, k consists of an absorption coefficient ic, and a 

6 o 

scattering coefficient 

"e ■ "a * 's • 

As a first order estimate of usually is ignored because its 

computation involves the dielectric properties of the ice particles and their 
size and shape distributions relative to the wavelength. For wet snow. 
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setting k- » Isa fair approximation because ic. « k. under most snow 
conditions In the microwave region. The $an» approximation holds true for 
dry snow at frequencies In the 1(M«er part of the microwave spectrum 
(< 5 GHz), but for frequencies above 10 GHz, may be comparable In 
magnitude to k, for some dry snow conditions (large particle sizes). Thus, 

O 

when the penetration depth Is calculated from &_ - 1/k., the values then 

P a 

obtained should be regarded as upper limits since they do not account for 
scattering losses In the medium. 

The absorption coefficient 1s defined as: 


tc 


a 




Xo 


- t(i + (4)')* 

2 e' 


1 ] 




where X^ Is the wavelength In free space. If fr-« 1. the above expression 
reduces to 


1C S — ■ ■ 

which Is applicable to dry snow throughout the micrwave region and 

to wet snow for V . < 2 %. 

w 

Figure 9 shows computed values of 6p plotted as a function of for 

three microwave frequencies. The values corresponding to dry snow 

conditions (V^ = 0) were computed using = 1.71 (frcxn Equation 6) and 

e"dj from Equation 7 with = 0.4 g cm"*and e"^ as given in Figure 3. 

For > 0, Equations 9-11 due to Linlor [31] were used for computing 

c'we ^"wc* l^^'ese equations were used at all three frequencies even 

though they are based on measurements over the 4-12 GHz region only. 

Furthermore, as noted earlier, because Linlor's measurements are representative 

of spherical Ice particles with coated layers of water which exhibit lower 

values of c".„ than snow under natural conditions, the values of 6^ shown 
ds p 

In Figure 9 are larger than would be expected for natural snow. That Is, 
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these curves should be regarded as estimates of the upper limit of 

Perhaps the most significant Information one can deduct from 
Figure 9 Is that the penetration depth can decrease by more than an 
order of magnitude as Increases from zero to 0.5X, followed by a much 
slower rate of decrease as a function of Increasing 

6.0 ENGINEERING APPLICATIONS 

Two of the uses of dielectric properties of snow are: (1) development 
of sensors to measure snow liquid water content and (2) Interpretation and 
modeling of microwave remote sensing data. For liquid water determination, 
frequencies below 100 ffiz have been found useful, while for remote sensing 
observations only frequencies above 1 GHz are usually employed. 

Knowledge of the liquid water content of snow Is of Interest In 
many areas of snow research. Including snow strength, snow melt, rate of 
metamorphism, water flow In snow, avalanche prediction and as a ground- 
truth factor for remote sensing observations at the higher frequencies. 
Airtoach and Denoth [25] have worked for several years to develop a dielec- 
tric measuring device that can be used to ascertain the liquid water 
content of snow. This Indirect method has numerous advantages over 
traditional methods such as calorimetric and centrifugal separation, 
among which are rapidity of measurement and reduction In the required 
ancillary equipment (dry Ice, silicone, etc.). Also when evaluated 
against the freezing calorimeter as a reference, an absolute accuracy 
of ±0.5% for was observed [35], This accuracy Is comparable to that 
of a field calorimeter and, although contact with the snow Is required 
for measurement, the time savings represented by this In-situ method, as 
compared with calorimetric methods, are very Important in field operations. 
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This technique Is based upon the substantial difference In the mag- 
nitude of the dielectric properties of water and dry snow as discussed 

In the preceding section. Ambach and Denoth's device consists of a 
bridge circuit that measures the capacitance of a probe filled with a 

snow sample. The change In capacitance due to the presence of the snow 
sample Is used to compute c* of the snow sample. Using Equation 10, and 
either measuring or estimating the snow density p^» the snow liquid water 

content Is obtained. In one configuration, the probe consists of 
5 parallel plates with a 3 cm sparing and an effective volume of approxi- 
mately 1000 cm. This volume was found to be near the lower practical 
limit based on constraints of minimal compaction of the snow sample and 
negligible fringing fields. 

In many cases liquid water determination In thinner layers Is de- 
sired; near the surface of a snowpack, the liquid water content can 
have a large gradient with depth. Planar comb-shaped capacitors have 
been developed [25] that, through control of widths and spacings, have 
varying effective measurement depths. Liquid water contents in layer 
thicknesses as small as millimeters can conceptually be obtained. Hence, 
the potential exists for measuring a very detailed liquid water profile 
of the snowpack near-surface layer which would be very useful In the 
interpretation of microwave remote sensing observations. 

Microwave radiometers and radars have shown promise for mapping 
snow depth, or water equivalent, over large areas for hydrologic appli- 
cations [32-34]. With the Ice crystals In snow being comparable In size 
to the wavelength in the microwave region, particularly In the milli- 
meter wavelength range, propagation In snow Is characterized by both ab- 
sorption and scattering. Hence the radar backscatter and radiometric 
emission from snow are ^vemed by the mean dielectric constant of the 
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snow volume as well as the spatial variability of the dielectric constant 
within the volume. Dry snow (and wet snow at low microwave frequencies 
In the neighborhood of 1 GHz) Is a relatively low-loss medium, which 
means that the observed backscatter and emission frmi a snowpack usually 
consist of contributions from the entire pack as well as from the ground 
below* Thus detailed knowledge of the snow structure and of the dielec- 
tric depth profile Is the key to successful modeling of the microwave 
Interaction with snow and to accurate Interpretation of microwave obser- 
vations. Due to the lack of experimentally-established models that can 
adequately describe the dielectric behavior of snow as a function of 
frequency and liquid water content, most research activities related to 
micrtmave remote sensing of snow have been limited to "first order" In- 
vestigations. In these first order Investigations, the mean dielectric 
constant usually is assumed to be uniform with depth and the snow crystals 
are assumed spherical In shape and uniform in size. These assumptions 
are used in theoretical and eir^lrlcal models developed to Interpret ex- 
perinental results. Such Investigations have established the overall 
behavior of backscattering and emission from snow, but In order to use 
remotely- sensed data to predict snow parameters— such as water equivalent 
and liquid water content— and to assess the state of the underlying 
ground medium (frozen or thawed), It will be necessary to use multi- 
frequency and multi -temporal observations In a prediction algorithm that 
Incorporates In greater detail the role of frequency and dielectric pro- 
file shape than Is currently possible. Thus, until further research of 
the microwave dielectric properties of snow Is conducted, development 
of remote sensing applications of snow can only proceed at a slow pace. 
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7.0 CONCLUDING REMARKS 


The electrical properties of snow have been evaluated as functions 
of the following ohyslcal properties: snow density, temperature. Ice 
particle shape, liquid water content, and the shape of the liquid water 
Inclusions. Analyses of the dielectric dependence on the above parameters 
Indicates that the liquid water content exercises the single roost Important 
influence, and consequently the discussions were divided Into wet and dry 
snow sections. Additionally, two spectral regions were examined In terms of 
potential applications. 

In the Irwer of the two regions, 100 KHz to 1 GHz, dielectric measure- 
ments have been used to measure the liquid water content of snow. By 
further restricting this range to be above 10 MHz, the dependence of 
on Ice particle shape becomes negligible. Hence, for a specific snow sample, 
only snow density must be measured for a calculation of liquid water content. 
Although not available commercially, operational sensors have been built, 
are in use for measuring liquid water in the snow, and hopefully will provide 
a good alternative to traditional techniques. 

In the spectral region applicable to remote sensing, 1 GHz to 35 GHz, 
the general electrical characteristics o'" snow are known; however addi- 
tional measurements are needed especially for the imaginary part of the 
dielectric constant above 10 GHz. Application of this information Is not 
direct as in the liquid water measurement, but it Is needed for refine- 
ment of the models used to describe the electromagnetic wave Interaction 
with snow. The sensitivity of microwave sensors to changes In the 
dielectric properties of snow (due to changes In liquid water content) may 
be illustrated by the following example. Corresponding to a change 
in liquid water content of the snow surface layer from = 0 (dry snow) 
to = 2%, the radar backscattering coefficient has been observed to change 
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by a factor of 30 (15 dB) at 35 GHz and the radlanetric temperature 
Increased by 120 K [35]. These results indicate that dry/wet snow 
determination is a relatively simple task, but quantifying the wetness ar 
measuring the water equivalent of the snowpack is not as straightforward 
To realise these latter objectives with a usable degree of prediction accuracy, 
improved knowledge of the microwave dielectric properties of snow is needed. 
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TABLE 1 


Measured Values of the Real Part of the Relative 
Dielectric Constant of Pure or Fresh-Water Ice, 

(from Ulaby et al. Lsi»}) 


Frequency 

(GHz) 

Temperature 

Range 

(“0 

^'l 

Reference 

0.15 to 2.5 

-1 to -60 

2.90 to 2.95 

Westphal (In [1]) 

9.375 

0 to -18 

3.15 

Gumming [14] 

10 

-12 

3.17 

Von Hippel [20] 

10 

0 to -35 

3.14 

Vant el al. [16] 

10 

-1 to -49 

3.17 

Lamb [17] 

24 

0 to -185 

3.18 

Lamb and Turney [18] 

26.4 to 40.0 

0 to -35 

2.92 

Vant et al. [16] 

94.5 

-28 

3.08 

Perry and Straiton [15] 

1000 

-173 

3.20 

Bertie et al. [21] 

99 

— 

3.17 

Blue [19] 
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Figure 2 , Cole-Cole diagram of e" and e* for ice samples at -10.8°C, 
a) pure ice and b) ice with impurities. Points represent 
measurements at the indicated frequencies in KHz. (From 
Evans [1]). 
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Figure 4. Dependence of the real part of the complex dielectric constant and loss tangent 
of dry snow with frequency for a) fine-grained and b) coarse-grained samples 
(from Ambach and Denoth [25]). 


Formula Frequency 

(Howorka lln31 1+2.22 P, 3.5 MHz 

\ Ambach And Denoth (251 1+2.205 P. 13.6 MHz 

— Kurloual29UEstlmated) 1+2.3 P, 3 MHz 

— Looyenga Formula II n 31 (1+0.508 P,) 



Figure 5. Real part of the dielectric constant of dry snow as a 
function of snow density along with empirical fits. 
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Figure 6. The loss tangent of dry snow as a function of temperature 
at a frequency of 9,375 GSiz (from Gumming [14]). 
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Pendular Regime— H \*- Funicular Regime 
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Figure 7. Depolarization factor A , for mixing Fomula 7 (from Ambach and Oenoth [25]). 
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Frequency; 9.4 GHz 
Density: 0. 5 g/cm^ 

- Measured Tobarias et al. (301 « Teniperature; 0**C 

- Calculations From Formula 6 
Calculations From Equation 10 &11 
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Figure 8b. Real and imaginary parts of the relative dielectric constant of snow as a function 
of liquid water content at 9.4 GHz. 









